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We review recent NLO QCD results for W, Z + 3-jet production at hadron colliders, computed using BlackHat 
and SHERPA. We also include some new results for Z + 3-jet production at the LHC at 7 TeV. We report new 
progress towards the NLO cross section for W + 4-jet production. In particular, we show that the virtual 
matrix elements produced by BlackHat are numerically stable. We also show that with an improved integrator 
and tree-level matrix elements from BlackHat, SHERPA produces well-behaved real-emission contributions. 
As an illustration, we present the real-emission contributions — including dipole-subtraction terms — to the pr 
distribution of the fourth jet, for a single subprocess with the maximum number of gluons. 



1. Introduction 

In the coming years a major theoretical task 
will be to provide reliable predictions for hard- 
scattering processes at the LHC. The start of the 
LHC era in particle physics opens new opportu- 
nities to confront data with theoretical predic- 
tions at scales well beyond those probed in pre- 
vious colliders. The first observation at the LHC 
of vector-boson production marks an important 
milestone. Processes involving vector bosons in 
association with jets are central to the physics 
program of the LHC. They are backgrounds to 
Higgs and top physics, as well as to many signals 
of new physics. Theoretical predictions can play 
an important role in experimentally-driven deter- 
minations of backgrounds, improving extrapola- 
tions to signal regions. Z and W production can 
play complementary roles in such determinations 
of backgrounds. At low luminosity, Reproduction 
can be used to calibrate estimates of Z produc- 
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tion, because its cross section is higher; while at 
high luminosity, Z production can be used to cal- 
ibrate W production because of the cleanliness 
of measuring lepton pairs. The good theoretical 
understanding of vector-boson production also al- 
lows its use to measure the luminosity. 

Next-to-leading-order (NLO) QCD is a key tool 
for confronting theory with experiment. This or- 
der in perturbation theory is the first to provide 
quantitatively reliable results, as it is the first 
order at which quantum corrections compensate 
the renormalization- and factorization-scale de- 
pendence in the strong coupling ofg. That scale 
dependence in leading-order (LO) predictions in- 
creases with the number of jets because of the 
increasing number of powers of as- At LO the 
scale dependence can be quite large, on the order 
of 50% for three or more jets. Furthermore, LO 
results may not model the shapes of distributions 
correctly. Although NLO computations are more 
challenging, in general they yield results with bet- 
ter reliability and agreement with measurements 
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(see e.g. refs. |ll2l3l4l5l6l7l8l^ l 

For many years the bottleneck blocking NLO 
computations of vector-boson production with 
three or more jets has been the difficulty of eval- 
uating the required one-loop amplitudes. The 
unitarity method [lOlllj . along with various im- 
portant developments j l2|13j . has broken past 
this barrier, producing NLO results for vector- 
boson production in association with three jets O 
[S]. (See refs. [416] for other unitarity-based re- 
sults, using various leading-color approximations 
to W+ 3 jets at NLO.) NLO computations of six- 
point processes involving heavy quarks have also 
been carried out recently, using both Feynman 
diagrams and on-shell methods |1]. 

In this talk we briefly review the NLO results 
for W, Z-|-3-jet production presented in refs. |5|8j . 
We also present new results on Z boson produc- 
tion at the LHC, all at the current center-of-mass 
energy of 7 TeV. Finally, we demonstrate the fea- 
sibility of computing W + 4-jet production at 
NLO. The latter process is key for new physics 
searches at the LHC, and from a theoretical point 
of view goes one jet beyond the 2007 Les Houches 
"experimenters' wishlist" [T]. We show that the 
virtual matrix elements computed by the Black- 
Hat library are numerically stable for the most 
complicated partonic subprocesses, those involv- 
ing the maximum number of gluons. We also 
show that, after augmenting SHERPA with tree- 
level matrix elements from BlackHat and in- 
corporating an updated integrator [M], we can 
obtain smooth distributions for the real-emission 
terms. As an illustration we present the infrared- 
subtracted real-emission contribution for a single 
subprocess. 

2. Calculational Setup 

NLO corrections to cross sections require the 
evaluation of virtual and real-emission contribu- 
tions. We evaluate the virtual amplitudes us- 
ing the BlackHat library |15I5) . a numerical 
implementation of on-shell methods |10I11I12I13| 
at one loop. For the real emission contribu- 
tions we use the Catani-Seymour dipole subtrac- 
tion method |16] , as implemented [17J in the pro- 
gram AMEGIC++ [18; (part of the SHERPA 



framework [TO]) to cancel the infrared divergences 
arising in the phase-space integration. For sub- 
processes with up to eight external legs we use 
AMEGICH — h to generate any required tree am- 
plitudes. We also use the SHERPA package to 
carry out the numerical integration over phase 
space. Details of the setup can be found in 
refs. |15I5I8J . 

3. NLO Results for Vector-Boson Produc- 
tion with Three Jets 




Figure 1. The scale dependence of the LO (dashed 
blue) and NLO (solid black) cross sections for Z, 7* -|- 
1,2,3-jet production at the LHC running at 7 TeV, 
as a function of the common renormalization and fac- 
torization scale fi, with fio = Mz. The bottom panel 
shows the K factors, or ratio between the NLO and 
LO result, for each of the three cases: 1 jet (dot- 
dashed red), 2 jets (dashed blue), and 3 jets (solid 
black) . 



In ref . ^ we presented a detailed study of !¥ -I- 
3-jet production at NLO, at the Tevatron and 
LHC. This study exhibited the expected scale- 
dependence reduction, and correspondingly bet- 
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ter theoretical precision, compared with LO pre- 
dictions. In addition, we found that at the LHC, 
the W~^ and W~ are both preferentiaUy polar- 
ized left-handed at large transverse momentum. 
This effect is also present in W^ -|- 1, 2-jet produc- 
tion, and at LO as well (see also ref. [20]). Such 
polarization will be largely, if not completely, ab- 
sent in W bosons emerging from decays of top 
quarks or other new heavy states. Accordingly, 
it should be useful in setting experimental cuts 
that will distinguish between such daughter Ws 
and "prompt" Ws emitted directly in the short- 
distance process. In this Proceeding, we will dis- 
cuss the more recent NLO computation oi Z + 3- 
jet production at the Tevatron [5], for which data 
already exists |3I7| . as well as new NLO predic- 
tions for the current LHC run. 

The dependence of NLO cross sections on the 
common renormalization and factorization scale, 
A* = f^R = Mf, is greatly reduced with respect to 
LO, as illustrated in fig. [1] for Z,j* + 1,2,3-jet 
production at the LHC at 7 TeV. For these plots 
the SISCONE jet algorithm [21] is used with cone 
size R — OA and merging parameter / — 0.75; 
the lepton and jet cuts are indicated on the fig- 
ure. (The case of three jets uses a leading-color 
approximation with the same separation of lead- 
ing and subleading terms described in ref. [8] ; we 
expect it to be accurate to within a few percent.) 
Indeed, the scale-dependence reduction becomes 
more pronounced as the number of jets increases, 
in accord with the increasing variation at LO due 
to increasing powers of as(^). The corresponding 
plot for production at the Tevatron shows similar 
behavior ^. 

In the left panel of fig. [5] we compare the dis- 
tribution in pt for the third-hardest jet {px- 
ordered) in Z, 7* -t- 3-jet production, computed 
at LO and NLO [8], to DO data. DO provided 
results for two data selections [7], with the Z 
decaying into and electron-positron pair. Fig. [5] 
uses the primary data selection based on the cuts, 
pi^* > 20 GeV, |77J^*| < 2.5, and 65 GeV < M^e < 
115 GeV. DO defined jets using the DO Run 
II midpoint jet algorithm [22], with a cone size 
oi R = 0.5 and a merging/splitting fraction of 
/ — 0.5. Because this algorithm is not infrared- 
safe, we use instead the SISCone algorithm [21], 



with the same parameters, and a scale choice 
fi ~ Ht/2. We also included non-perturbative 
corrections estimated by DO [7]. 

The figure demonstrates good agreement be- 
tween NLO theory and experiment, despite a 
number of issues, including the difficulty of ob- 
taining precise nonperturbative corrections, and 
differences in the jet algorithms. Further details 
of comparisons to CDF and DO data for Z-boson 
production in association with jets may be found 
in ref. 8J. 

The right panel of fig. [2] shows the same distri- 
bution at the LHC (except that it is not normal- 
ized to (7z,'y')- The experimental cuts are labeled 
on the plot. We again use the SISCone jet al- 
gorithm, which (along with anti-Zcy) is one of the 
infrared-safe jet algorithms adopted by the LHC 
experiments. Our setup allows results for differ- 
ent algorithms and cone sizes to be computed at 
the same time. It should be possible to minimize 
non-perturbative effects by choosing appropriate 
cone sizes [23] , perhaps in a pT-dependent fash- 
ion. 

4. Towards W + A Jets at NLO 

W + 4-jet production has been an important 
background since the early days of the Tevatron. 
There, it was the dominant background to ti pro- 
duction; at the LHC, it will be an important back- 
ground to many new-physics searches, and will 
continue to be important to precision top-quark 
measurements. Its computation beyond LO rep- 
resents an important challenge to theorists. We 
report here on progress towards an NLO calcula- 
tion of this process. 

A crucial issue in the approach pursued in 
BlackHat is numerical stability of the virtual 
contributions. In fig. [3] we illustrate the stabil- 
ity of the leading-color virtual corrections to the 
squared matrix element, day, summed over colors 
and over all helicity configurations for the subpro- 
cess gd -^ e~vgggu. The horizontal axis of fig. [3] 
shows the logarithmic error. 
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for each of the three components: 1/e^, 1/e and 
e'^, where e = (4 — D)/2 is the dimensional reg- 
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Figure 2. The left plot shows l/az,-y» x da/dpT for the third-hardest jet in Z, 7* + 3-jet production, comparing DO 
data against LO and NLO predictions. In the upper panel the parton-level NLO distributions are the solid (black) 
histograms, while the NLO distributions corrected to hadron level are given by dash-dot (magenta) histograms. 
The DO data is indicated by the (red) points. The LO predictions corrected to hadron level are shown as dashed 
(blue) lines. The lower panel shows the distribution normalized to the full hadron-level NLO prediction. The 
scale-dependence bands in the lower panels are shaded (gray) for NLO and cross-hatched (brown) for LO. The 
right plot shows the pr distribution for the third jet at the LHC, with the indicated cuts. The complete color 
dependence is included in this plot, but no corrections to hadron level are applied. 



ularization parameter. In this expression ay is 
the cross section computed by BlackHat as it 
normally operates for production runs, whereas 
^target .^ ^ target value computed by BlackHat 
using higher-precision arithmetic of at least 32 
digits. The phase-space points are selected in the 
same way as those used to compute cross sec- 
tions. As seen in fig. El no point has a relative 
error larger than about 1%; this level of precision 
is more than adequate for phenomenology. 

What about the real corrections? In fig. |4] we 
show the real-emission contributions, including 
dipole subtractions, for the distribution in the 
fourth-jet pr- Although these contributions alone 
are not physically meaningful (they need to be 
combined with the analytically-integrated dipole 
terms and the virtual corrections, and at the very 
least summed over other single-quark-line sub- 
processes), they do serve to illustrate that our 
SHERPA-based integration setup has them un- 
der control. Here we used the anti-fc^ algorithm 



for R ~ 0.4, with a jet cut oipT > 25 GeV, along 
with other cuts on the leptons and the jet rapidi- 
ties. For this plot, the real terms were sampled 
at 2 X 10* phase-space points, yielding the rather 
small integration errors indicated by the thin ver- 
tical lines. Two key improvements in our setup, 
which make the real contributions feasible, are an 
improved integrator [1418] and more efficient nine- 
point tree-level amplitudes generated by Black- 
Hat using on-shell recursion relations [24) . 

5. Conclusion 

A publicly available version of BlackHat is 
under construction and is being tested in a va- 
riety of projects (see e.g. ref. [25]). This version 
uses the Binoth Les Houches interface [26] . It has 
been tested with both C++ and Fortran clients. 
We intend the public version to provide all pro- 
cesses that have been thoroughly vetted using the 
full BlackHat code. We expect significant fur- 
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Figure 3. The distribution of the relative error in the 
leading-color virtual cross section for the subprocesses 
gd -^ e~i'gggu. The horizontal axis is the logarithm 
of the relative error between an evaluation by Black- 
Hat, running in production mode, and a higher pre- 
cision target. The vertical axis shows the number of 
phase-space points out of 100,000 that have the cor- 
responding error. The dashed (black) line shows the 
1/e^ term; the solid (red) curve, the 1/e term; and 
the shaded (blue) curve, the finite (e°) term. 
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Figure 4. The real part of the NLO pr distribu- 
tion for the fourth jet in VK -I- 4-jet production at 
the LHC. Only the subprocess gg -^ e~DggQqg is 
included here. The thin vertical lines, where visi- 
ble, indicate numerical integration uncertainties. 



ther gains in efhciency after implementing a va- 
riety of additional optimizations, such as taking 
advantage of phase-space symmetries to reduce 
the amount of computation. 

In summary, V + 3-jet production (where V 
is a vector boson) at hadron colliders is under 
excellent control at NLO, using BlackHat and 
SHERPA, allowing comparisons [4|5|8j to Teva- 
tron data [21317) and predictions for the LHC. The 
frontier has now shifted to ^ -I- 4-jet production. 
In this Proceeding, we presented significant new 
progress towards this goal, and demonstrated for 
a selected subprocess that the programs are able 
to execute these complex calculations. After fur- 
ther optimization and testing we plan to compute 
V + 4-jet production at NLO at the LHC. More 
generally, we look forward to applying these tools 
to a wide range of studies of LHC physics. 
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